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1 I The acoustic properties of a porous sheet of medium resistivity backed by a rigid 

plate in which are embedded a periodic set of circular inclusions is investigated. 

O . Such a structure behaves like a multi-component diffraction gratings. Numerical 

O ' results show that this structure presents a quasi-total absorption peak below the 

• i-H ' 

£>y essential spectrum, i.e. below the frequence of the fundamental quarter-wavelength 

Q_i! resonance of the porous sheet in absence of inclusions. This result is explained either 



by the excitation of a complex trapped mode, or by the increase of viscous loss 
associated with a larger velocity gradient inside the layer at the modified quarter- 
wavelength resonance frequency. When more than one inclusion per spatial period 
are considered, additional quasi-total frequency peak are observed. The numerical 



l> 

results as calculated with the help of the mode-matching method described in the 



paper agree those as calculated with the help of a Finite Element method. 
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I. INTRODUCTION 

This work was initially motivated by a design problem connected to the determination of 
the optimal profile of discontinuous spatial distribution of porous materials and geometric 
properties for the absorption of sound. Porous materials (foams) suffer from a lack of 
absorption at low frequencies, when compared to their efficiency at higher frequencies. The 
purpose of the present article is to inverstigate an alternative to multi-layering by considering 
periodic inclusions embedded in a porous sheet attached on a rigid plate. This configuration 
results in a diffraction grating and possibly a sonic-crystal used in reflection. 

The influence of a volumic heterogeneity addition on absorption and transmission of a 
porous layer was previously investigated by use of the multipole method ir*i>2, by embedding 
a periodic set of high-contrast inclusions, whose sizes are comparable to the wavelength, 
in a macroscopically-homogeneous porous layer whose thickness and weight are relatively 
small. This leads either to an increase of the absorption coefficient in the case of one 
layer of inclusions or to band-gaps and a total absorption peak in case of multi-layered 
set of inclusions (sonic crystal). The influence on the absorption was explained by mode 
excitation of the configuration enabled by the periodic inclusions, whose structure leads 
to energy entrapment. Other works related to volumic heterogeneities in macroscopically 
homogeneous porous material were carried out essentially by means of the homogeneization 
procedure^, and possibly leading to double porosity materials^. 

The influence of the irregularities of the rigid plate on which are often attached porous 
sheets on the absorption coefficient was previsouly investigated by use of the multi-modal 
method in^, by considering periodic rectangular irregularities filled with air. This partic- 
ularly leads in the case of one irregularity per spatial period to a total absorption peak 
associated with excitation of the fundamental modified mode of the backed plate. The lat- 
ter is excited thanks to the surface grating. Other works related to surface irregularities 
were carried out, notably related to local resonances associated with fractal irregularities^. 

Local resonance and trapped modes are an other possibility to localize the field. Trapped 
mode were largely studied in guided waves^ or in periodic structures^. Here, we investigate 
theoretically and numerically the influence on the absorption coefficient of the embedding 
of multiple inclusions grating in a rigid frame porous layer glued against a rigid wall. The 
effects of the modified mode of the plate and Bragg interference are clearly visible on the 




FIG. 1. Cross-sectional plane view of the configuration. 

absorption curve, while a quasi-total absorption is obtained for a frequency below the one of 
the fundamental quarter- wavelength resonance of the backed porous sheet. The latter peak 
presents some of the specific features of a trapped mode excitation. 

II. FORMULATION OF THE PROBLEM 
A. Description of the configuration 



Both incident plane acoustic wave and plate are assumed to be invariant with respect 
to the Cartesian coordinate X3. A cross-sectional x\ — 22 plane view of the 2D scattering 
problem is shown in figure [TJ 

Before the addition of the cylindric inclusions, the layer is made of a porous material 
saturated by air (e.g., a foam) which is modeled (by homogenization) as an (macroscopically 
homogeneous) equivalent fluid M'- 1 '. The porous sheet is backed by a rigid surface. The 
upper and lower flat and mutually-parallel boundaries of the layer, whose 22 coordinates 
are H and 0, are designated by Th and Tq respectively. M' ' and ilfM are in firm contact 
through T H , i.e. the pressure and normal velocity are continuous across T H ([p(x)] = 
and [p~ 1 d n p(x.)] = 0, wherein n denotes the generic unit vector normal to a boundary and 
d n designates the operator d n — n • V). T is rigid (Neumann type boundary conditions, 
<9 n p(x) = 0). 

N c inclusions with a common spatial periodicity d are embedded in the porous layer, 
that create a diffraction grating in the x\ direction. Depending on the arrangement of 
the N c inclusions in the unit cell, a diffraction grating or a sonic-crystal of period d c can be 
formed (d c < d). The set of indices by which the cylinders within the unit cell are identified is 



denoted by J\f c E N. The j'-th inclusion occupies the disk f2' 2 3 ' of radius R^ and is centered 
at x^) = (xi , X2 ), J S A/" c . The inclusions are infinitely rigid (Neumann type boundary 
conditions on r^'), i.e. the contrast between the elastic material M™ and M^ 1 ' is very large. 
This also means that the inclusions can consist in tubes or holes posteriorly proofed for 
acoustic waves. Two subspaces fi' 1 ' G f^ 1 ' are also defined respectively corresponding to 
the upper and lower part of the plate not without inclusions. 

The total pressure, wavenumber and wave speed are denoted by the generic symbols p, 
k and c respectively, with p = p[°\ k = k^ = u/c^ in $7^ and p = p^\ k = fcM = ui/c^ in 

Rather than to solve directly for the pressure p(x, t) (with x = (xi,X2)), we prefer to 

/OO 
p(x,cu)e~ luit duj. 
-oo 

Henceforth, we drop the u in p(x, u) so as to denote the latter by p(x). 

The wavevector k 1 of the incident plane wave lies in the sagittal plane and the angle 
of incidence is 9 l measured counterclockwise from the positive x\ axis. The incident wave 
propagates initially in f^°l and is expressed by p*(x) = A t e 1 ^ xi ~ k2 ( :E2_ - ff )) ) wherein k\ = 
— fc[°] cosO 1 , k 2 = fc'°l sin#* and A 1 = A l (uS) is the signal spectrum. 

The plane wave nature of the incident wave and the periodic nature of II £T 2 3 ' imply 

jeN c 
the Floquet relation 

p( Xl + qd, x 2 ) = p(x 1 ,x 2 )e ik ^ d ;VxGl 2 ;VgGZ. (1) 

Consequently, it suffices to examine the field in the central cell of the plate. 
The uniqueness of the solution to the forward-scattering problem is ensured by the radi- 
ation condition: 

p'°l(x) — p % (x) ~ outgoing waves; |x| — » oo, x 2 > H . (2) 



B. Material modeling 

Rigid frame porous material M is modeled using the Johnson-Champoux-Allard model. 



The compressibilty K and density p, linked to the sound speed through c = y/TJ (Kp) 



aj.gjjj.il 



K i / / \ - 1 



1 T^b 

7-(7-l)fl + i^G(Pr^ ^ (3) 

wherein cu c = acj)/ Pfa^ is the Biot frequency, u' c = a'tfi/pfOtoo, 7 the specific heat ratio, 
Pq the atmospheric pressure, Pr the Prandtl number, pj the density of the fluid in the 
(interconnected) pores, the porosity, a^ the tortuosity, a the flow resistivity, and a 1 the 
thermal resistivity. The correction functions G(Pra;)^ and F(oj)^ are given by 



0(P»)=^/l-^/P»(^) , 



F(w) = Wl — ir/pfiu 



o \ 2 



(T0A / 

where 77 is the viscosity of the fluid, A' the thermal characteristic length, and A the viscous 
characteristic length. The "thermal resistivity" is related to the thermal characteristic length 
through a' = 80^1] /cp A' 2 ,—. 

The configuration is similar to but also more complex than those already studied in^»^, 
in sense the porous sheet is backed by a rigid plate, but also in sense the unit cell can be 
composed of more than one non overlapping inclusion, when the a^-coordinates of the center 
of two inclusions are separated by a distance lower than the sum of their radii. In this latter 
case, the interaction between these inclusions can not be modeled as exposed in^ and a more 
complex interaction model should be employed^. The method of solution is also briefly 
summarized hereafter. 

C. Field representations in fij® and fi^ 1 1 

The continuity relations across the interfaces Th and To are first considered in section 
IIII Al The field representations in DJ-°\ and f^ 1 1 are needed as the first step. The continuity 
conditions across T^\ \/j E J\f c will be treated in section UlIBI 

Separation of variables, the radiation condition, and the Floquet theorem lead to the 
representation: 

p[°](x) = J2 [e- ik ™ (x2 - H) 5 q + R q e<^- H A e [k ^\ Vx e ^ , (5) 



wherein 5 q is the Kronecker symbol, k\ q = k\ + 2qir/d, k 2q = \j (k^) 2 — (ki q ) 2 , with 
Re ( k 2q ) > and Im (k 2q ) > 0. R q is the reflection coefficient of the plane wave denoted 
by the subscript q. 

It is first convenient to use Cartesian coordinates (xi,X2) to represent the field in 
fit 1 1. The latter are composed of the diffracted field in the plate and the fields scat- 
tered by the inclusions, whose form depends on the position of x, either below or above the 
inclusions 16 . Refering to^, whatever the arrangement of the inclusions, £2 is always larger 
than maxjg^c (x 2 + Ry> ) in fit 1 ', while 22 is always smaller than imuj^j^c (x 2 — R^' ) 
in fit 1 1. The total field in fit 1 1 can be written in Cartesian coordinates as 



J2(f q e-<^+ 9q e 



p* f (x)= > lt q e"**"* + g q (r**-*)e 



E E E tf^VM 21 ^^ h-4 j) )) , ( 6 ) 



q&j£Af c l£2 



wherein B\ are the coefficients of the field scattered by the j'-th cylinder of the unit cell, f q 



and g q are the coefficients of the diffracted waves inside the layer associated with the plane 

:± = 2(-i) l e ±iw "/dk [ 2q ] with 6 q such that kV« = k lq + &£]. 



wave denoted by q, and K± = 2(-i) 1 e ±ue " /dk [ 2q ] with 6 q such that k^e ie " = k lq + ik l 2q \ 15 . 



III. DETERMINATION OF THE ACOUSTIC PROPERTIES OF THE 
CONFIGURATION 

A. Application of the continuity conditions across Th and To 

Applying the continuity of the pressure field and of the normal component of the velocity 
across Th and the Neumann condition on To, introducing the proper field representation 

therein, Eqs.(j5]) and ([6]), and making use of the orthogonality relation / e 1 ^ ln ~ u ' xi dxi = 

J-l 

d5 n i, V(/,n) e 1? give rise to a linear set of equations. After some algebra and rearrange- 
ments, this linear set reduces to a coupled system of equations for solution of R q , f q and g q 
in terms of B\ . 



B. Application of the multipole method 



The expressions of f q and g q in terms of B) 3 ' are introduced in the so-denoted diffracted 
field inside the layer. The latter field accounts for the direct diffracted waves inside the 
layer and for the reflected waves at the boundaries T and Th previously scattered by each 
inclusions. This expression when compared with the expression of the direct scattered field 
by the inclusions, is valid in the whole domain f)W. To proceed further, the Cartesian form 
of this field is converted to the cylindrical harmonic form in the polar coordinate system 
attached to each inclusion, as stated for example ir*i. Effectively, central to the multipole 
method are the local field expansion or multipole expansions around each inclusion. 

The pressure field in the vicinity of the J-th inclusion in the polar coordinate system 
attached to this inclusion and introducing A L the coefficient of the locally-incident field, 
reads as 



pW(r J ) = X)4 ,7) Hi 1) (*Wr J ) e 
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E ff«s. w + E E s^sf + e E E <W + E *; 
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qL 



lei 



j^J lei 



jeM c lei qei 



Jl (k [1] rj) P*' 



Y^lB^ik^+A^J^rj) 



Lei 



J-LBj 



(7) 



r qL 



O iJ ' j) 

^Llq 



2S ^ cos (k^-L9 q )e<^ 



D„ 



Sl-i= £h« (JfcW«J) 



i=l 



lk{id 



+ (-1) 



2(-iy- L e i k iM J) - x i y ) 

h 2p ±J P 



L-l -\k\id 



Ji]; 



dk^D. 

ii. ^mV,." 1 ,'■" 



(af 1 - 4°]) e^ H cos 4 1 ' (xf - x^ ) - (/ - L)9 q 



+af cos f fcg I x K i> + x ( 2 J> -H)-(l + L)6 q ) + ia[ 0] sin f fcjg f x^' } + x[ J) - H) - (I + L)6 q 

D q = af cos (kf q H\ - ia£ ] sin fe ] #) , 

(8) 



wherein Sl_/ is the lattice sum often refered to as Schlomilch series for non-dissipative 
material, H^ is the L-th order Hankel function of first kind and J^ is the L-th order Bessel 



function. The terms S L ^ accounts for the coupling between the multiple inclusions inside 
the unit cell and take the form: 

S]H> = £ (-i)"2e^-'>» V(t „, J . J ..,. i ., ± 4;i ( 4,..,o. )) (1 _ hi) (9) 

wherein the signs + and — correspond to x J 2 > x\ and x 3 2 < x\ respectively, which can be 
found in 1 when \xi — xi \ > R^ + R® or 



s Ll 



(Jj) 



BiMl (r}) e^^+^S^i^^^VJ e 1( °-^ 



>JiJ 



(10) 



when |x 2 — a?2 | > -R^ + -ft • This latter form agrees with the one found in^ when the 
inclusions are aligned inside the unit cell, i.e. x$ = x 2 , Vj G A^ c , which imposes # j = or 
#j = n. In Eg. (II Op. (rj, 9j) is the coordinate of x^ in the polar coordinate system attached 
to the J-th inclusions, i.e. centered at x( J ). 

Finally, it is well known that the coefficients of the scattered field and those of the locally- 
incident field are linked by a matrix relation derived from the boundary condition on r( J ) 
only, i.e., B L = V£ A L , wherein V£ are the cylindrical harmonic reflection coefficients. 
These coefficients takes the form V£ = — H L (k^'R^') /Jt, (k^'R^)) in case of Neumann 
type boundary condition, with x{ x ) = dx/dx, x being either Hankel or Bessel functions. 
Introducing the expression of A L derived from Eq. (|7|) in the previous relation gives rise to 
the linear system of equations for the solution of B[ . This linear system may be written 
in the matrix form, when denoting by B the infinite column matrix of components B[ 

(I-V(S + Q))B = VF, (11) 

wherein F is a vector of components y, FqL ■> which accounts for the solicitation of the J-th 

inclusion by a wave that is previously diffracted inside the layer, V is a diagonal square 
matrix of components V£ , S and Q are two matrices of components Sl-i + S^/, which 
accounts for the coupling between the J-th and the j-th inclusion inside the layer and 
/ jQliq ii which accounts for the coupling between the J-th and the j-th inclusion through 

waves diffracted by the layer. 

The expressions of the components involved in (I — V (S + Q)) B = VF are identical to 
those found in 1 , when the half-space behind the layer vanishes and when the center of the 
cylinders are defined as they are in the present article. 
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C. Evaluation of the fields, reflection and absorption coefficients. 

Once the linear system (TTTT) is solved, the expression of R p in terms of B\ 3 reads as 
Hcoste'^H-iaW^infe 11 ^) 



a 1 
R„ = — 



D' 



A t -\ l M 

qez jeAf c i& dk 2q u v 

The first term corresponds to the reflection coefficient in terms of waves in absence of 
inclusion, i.e. for q = or for the incident plane component indexed by i, and the second 
terms accounts for the inclusions. 

Introduced in Eq.(j5]), the reflected field is expressed as a sum of the field in the absence 
of the inclusions with the field due to the inclusions. In case of an incident plane wave with 
spectrum A l (ou), the conservation of energy relation takes the form 

I = 11 + A, (13) 

with TZ and A the hemispherical reflection and the absorption coefficients. TZ is defined by 



Re(k [ °}) iin i|2 «+ jlM no ,|2 

V 2q J \\ K q\\ V^ 2 1 \\ K lW 



' X Z^ , [0]i 11^112 Z^ , [0]i \Ui\\2 ' ^^ 

ft 2 II II q=-q- K 2 II I' 



wherein g T are such that g T < rf/27r (/c^ ± k\) < q T + 1 and the expression of _R P are given 
equation (fl2l . .4 takes the form 

^ = „ rm,- „ (^ + As) > ( 15 ) 



d^ UJ, ||A<" 2 



whererin 



^ = ReT^I) j[ m Im (^') 2 ) l|p[1] (X) l|2d " ' (16) 

is the inner absorption of domain fiM. dw is the differential element of surface in the sagittal 
plane and 

• 4s=Re 15MS p,l, * (x) "° i ' vp,l,(x)d7 ' (i7) 

is the surface absorption related to interfaces Th- dj is the differential arc length in the 
cross-sectional plane, z/ i is the outward-pointing unit vector to the boundary T H , and p* is 
the complex conjuguate of p. 



TABLE I. Geometry of the considered configurations. 





N 


d (cm) 


H (cm) 


(x-L (cm),x 2 (cm)) 


flCj) (cm) 


CI 


1 


2 


2 


(1,1) 


0.75 


C2 


2 


2 


3.5 


(x^^na,!) 

(xi 2) ,4 2) ) = (0.5, 1 + \/3d/2) 


itf 1 ) = 0.75 
i?( 2 ) = 0.5 



.As accounts for the absorption induced by the viscous dissipation at the interfaces. 
Effectively, it is obvious from Eq. (IT7|) that As does not vanish because of non- vanishing 
Im (p^ 1 '), which is a consequence of the modeling of viscous dissipation phenomenon^. 

Because of the complicated shape of fiW anc j the non- vanishing term As, A will not be 
calculated by the expression given in ffTol) . but rather by A = 1 — 1Z. 

IV. NUMERICAL RESULTS, VALIDATION AND DISCUSSION 



The infinite sum X]«ez over the indices of the k± q depends on the frequency and on 
the period of the grating. An empirical truncation rule is employed, inspired from^^ 
and determined by performing a large number of numerical experiments Ylq=-o_ such that 
Qzp = int {(1/271 (3Re (fcM) ± fc})) + 10. In the latter equations, int (a) represents the integer 
part of a. 

The infinite sum Xlmez over the indices of the modal representation of the diffracted field 
by a cylinder is truncated 18 as Ylt=-u such that M = int [Re (4.05 x (k^R) Ti + k^R\) + 
10. 

Finally, the infinite sum (lattice sum) embedded in Sl-i in equations (jSJ) and ffTUl) Y^=i i s 
found to be slowly convergent, particularly in the absence of dissipation, and is found to be 
strongly dependent on the indice L — I. A large literature exists on this problem^^. Here, 
the fact that the medium M"W is dissipative greatly simplifies the evaluation of the Schlomilch 
series. The superscript I in S^Jj identifies the integer over which the sum is performed, i.e. 
^2 i=v This sum is carried out until the conditions Re ( [S^ i _ l — S^IA / S^A < 10~ 5 
and Im ((<Sf_+ 1} - sg,) /Sff,) < 10 ~ 5 are reached^. 

Numerical calculations have been performed for various geometrical parameters whose 
values are reported in Table [H and whithin the frequency range of audible sound, particularly 



10 



TABLE II. Parameters of the porous foam used in the article. 



tt'r 



A (pm) 



A' (pm) 



a (Ns.m 



v c (Hz) 



0.95 



1.42 



180 



360 



8900 



781 



at low frequencies. For all calculations, the ambient and saturating fluid is air (p^ = 
p f = 1.213 kg.m" 3 , cM = ^Po/pf, with P = 1.01325 x 10 5 Pa, 7 = 1.4, and 77 = 
1.839 x 10~ 5 kg.m" .s _1 ). One of the main constraints in designing acoustically absorbing 
materials are the size and weight of the configuration. A particular attention is paid on 
the dimension, i.e. thickness, and the frequencies of the absorption gain, which have to be 
commonly as small as possible. The absorption gain is defined by reference to the absorption 
of the same configuration without inclusion embedded. The initial configuration consists in a 
2 cm thick porous sheet of Fireflex (Recticel, Belgium) backed by a rigid plate. The material 
characteristics are reported in Table HT1 and were determined by use of traditional methods 10 . 
Circular cylinders of 7.5 mm radius are embedded with a spatial periodicity of 2 cm. 

A. One inclusion per spatial period 

We first consider only one inclusion centered in the unit cell, i.e. (x[ ,x 2 ) = (d/2, H/2) = 
(1 cm, 1 cm). The first two modified modes of the plate, which are excited because of the 
periodic arrangement of the inclusions, Appendix |Al stand around ^(1,1) ~ 14 kHz and 
^(2,1) ~ 16 kHz. The attenuation associated with both modes is relatively large, figur^3j 

Different type of waves corresponds to each kind of mode related to the grating, i.e. 
mode of the grating (MG) and modified mode of the backed layer (MMBL): evanescent 
waves in fiM (and also in fi^l) for the MG, and evanescent waves in f^ ' and propagative 
waves in fiW for the MMBL. In order to determine which type of mode is excited by the 
plane incident wave, we have plotted in Figure [2] the transfer function as calculated by 
TF[y) = p(x, w)/V°' l (x, uS) on r (x 2 = 0) at 1 cm from the center of the inclusion (between 
two inclusions), when excited at normal incidence. The tranfert function is separated on the 
different intervals corresponding to the different type of waves that are involved in the total 
pressure calculation: TF(y) is the total transfer function, TF\{y) is the contribution of the 
propagative waves in both fl^ and V^ l \ TF 2 (u) is the contribution of the evanescent waves 
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FIG. 2. Configuration CI - Transfer function, TF ( — ), on Tq at 1 cm from the center of the 
inclusion (between two inclusions), and its different contributions when the configuration is excited 
at normal incidence : ( ) TF\, (—■) TF2, and (••) TF3. 

in f^°l and propagative ones in fiW, and TF 3 (u) is the contribution of the evanescent waves 
in both fi^ and fiM. The transfer function possesses a large peak at w 15 kHz, around 
^(1,1) and ^(2,i)- This also proves that the MMBL are the most excited modes related to the 
grating, at least around these frequencies. The peak results from a continuous drop between 
evanescent waves in both material to evanescent waves in the air medium. This also means 
that this peak is neither a MMBL nor a MG, but result from a complex combination of these 
two types of mode, with a structure closer to the one of the MMBL. Because of this stucture, 
the energy is trapped in the layer leading to an increase in the absorption coefficient. The 
translation of the excitation of these modes in terms of absorption, i.e. the peak around 
17 kHz figure [3J is smaller than the one depicted in^, because the attenuation associated 
with the modes in the present configuration is larger and because the flow resistivity of the 
foam considered here is larger. The design of a struture composed of a layer with inclusions 
embedded is more based on compromises than with irregularities of the rigid backing 2 . These 
compromises relate the spatial periodicity, the radius of the inclusion or better the ratio R/d 
which can not be too small and which constrains the layer thickness and the properties of 
the latter. 

Figure [3] depicts the absorption coefficient as calculated for this geometry. This result 
was validated numerically by matching the absorption coefficient as calculated with the 
present method with the one as calculated with a Finite Element method. Quadratic finite 
element were used to approximate the pressure inside the unit cell, thereby leading to a 
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discretized problem of 2196 elements and 1238 nodes. The periodicity relation, i.e. the 
Floquet condition, were applied on both sides of the discretized domain, i.e. at each nodes 
of Xi-coordinate and d. For this periodicity relation to be correctly implemented these 
two sides were discretized with similar nodes, i.e. idential a^-coordinate. The results match 
well, thus validating the described method, figure EJ 

Because of the rigid backing, which acts as a perfect mirror, the response of the con- 
figuration possesses some particular features related to multi-layered grating. We also in- 
troduce d 2 = 2x 2 the distance between the center of the circular cylinder and the center 
of its image. Each grating is going to interfere with one another at the Bragg frequencies 
u) n \ = rm~Re (c^) jdi- In particular, the first Bragg frequency is v h ,^. « 6 kHz. The lat- 
ter frequency is largely employed to determine the central frequency of the band gaps for 
phononic crystals and corresponds to the maximum of reflected energy and so to a minimum 
of transmitted energy -band gaps- in case of phononic crystal. The absorption coefficient 
also presents a minimum at v\y 

A particular feature of the response of this configuration is that the absorption coefficient 
presents a peak close to unity at a low frequency v t below the so-called fundamental quarter- 
wave resonance frequency, i.e. below what can be associated with an essential spectrum. 

A sensitivity analysis, performed by varying one parameter while the other are kept 
constant at their value, shows that the radius of the inclusion has a large influence on v t and 
on the amplitude of the corresponding absorption peak. The radius R was varied from 1 mm 
to 9.5 mm. In terms of amplitude of the absorption peak, R = 75 mm is the optimal value, 
while v t decreases when R increases. In the opposite, d = 2 mm is the optimal value in terms 
of the amplitude of the peak, but v t increases when d increases from 1.75 cm to 3.75 cm. The 
spatial periodicity of the arrangement acts inversely on v t than it does on the frequencies 
of the modes closely related to d like the MMBL and the MG, Appendix |Aj In fact, the 
amplitude of the absorption peak increases with the filling ration R/d until being close to 
unity and drastically decreases after this value because the wave can no more propagate in 
the layer towards the rigid backing and is mainly reflected on the cirular grating, or because 
the density of inclusion becomes insignificant. In the same way, when the angle of incidence 
decreases [tt/2; vt/6], the v t increases in the opposite to the frequencies of the MMBL and of 
the MG. The amplitude of the peak is quite close to unity until 7r/6. For smaller values of 
incidence angle, the amplitude of this low frequency peak begins to decrease. 
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FIG. 3. Configuration CI - Absorption coefficient of a H = 2 cm thick porous sheet of Firefiex 

backed by a rigid plate ( ) without inclusion embedded and ( — ) with a R = 75 mm radius 

circular cylinder embedded per spatial period d = 2 cm. The Finite Element result is plotted with 
(o). The absolute value of the determinant of the propagation matrix D = 1 — V(S + Q) is plotted 
on top of the figure. 



When R = 75 mm and d = 20 mm, v t is all the smaller that the inclusion is distant from 
the rigid backing, i.e. that x 2 , the center of the inclusion, is large. The amplitude of the 
peak is close to unity whatever x 2 , in [0.8 cm; 1.2 cm] the range conditionned by the layer 
thickness and the inclusion radius. 

The features of this low frequencies absorption peak sounds like phenomena that are 
related to trapped modes in wave guides^ or embedded Rayleigh-Bloch waves^ 1 ^. These 
modes has finite energy and corresponds to a solution which decays down away from the 
perturbation. Figured shows a snapshot of the module of the pressure field at v t . The 
latter clearly exibits a maximum on the side of the rigid plate and a minimum on the side 
of To, which is typical of a trapped mode. Everything seems to happen like if a Dirichlet 
wave guide of thickness 2H presents, necessarily, symetric obstacles formed by the inclusion 
and its image. In our case, these trapped modes are complex, because the boundaries of 
the wave-guide are not Dirichlet conditions but continuity conditions, and because M^ is 
a dissipative medium. The determinant as calculated for the configuration CI, figure |3j 
present a minima at u t , which suggests that a complex (trapped) mode CTM stands at 
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FIG. 4. Configuration CI - Snapshot of the module of the pressure field inside the porous sheet at 
u t = 2674 Hz. 

this frequency. In the opposite, it is clear from figure^ that the peak around v t is mainly 
associated with propagative waves in both domains, a small part of it being associated 
with evanescent waves in the layer which entraps the energy. This phenomena was already 
encountered in^ and attributed to the periodicty of the configuration. An other explanation 
of the quasi-absorption peak is related to the modification of wave path and structure global 
properties inside the porous sheet. For a particular ratio R/d, the pressure gradient, which 
is a cause of viscous loss, is larger inbetween the inclusions and between the inclusions 
and the rigid backing. This loss phenomenon results from a continuous modification of the 
quarter-wavelength resonance when the inclusion radius increases. 

Other layer thicknesses were tested. It was found that for each layer thickness, in the 
suitable range for the application, and a centered (x[ ,x 2 ) = (d/2, H/2) inclusion, a couple 
(R^',d) exists for which a quasi-total absorption peak exists below the quarter- wavelength 
resonance frequency. 

Finally, because the parameters of a foam are often difficult to predict before its polymer- 
ization, a sensitivity analysis has been performed with regards to the acoustic and structural 
parameters of the porous sheet. The amplitude and frequency of the CTM is quasi indepen- 
dant from a variation of <f), A and A'. When a increases from 1.02 to 1.42, the sound speed 
in the material decreases and v t decreases, while the amplitude of the associated peak stands 
close to one. The resisgtivity a particularly influences the amplitude of the peak. When 
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it increases, the amplitude admits a maximum and the peak is wider, while v t increasse. 
The resistivity a ([3900 Ns.m" 4 ; 12900 Ns.m -4 ]) is the parameter that mostly influences the 
results and its value has to be close to the one used in the simulations, i.e. in our case, a 
value inbetween 7000 Ns.m -4 to 11000 Ns.m -4 is acceptable. 

B. Two or more inclusions per spatial period 

Various configurations were tested, involving two or more inclusions per spatial period. 
The frequency band investigated stands below the quarter-wavelength resonance frequency 
of the associated porous sheet or at least below the first Bragg frequency, i.e. below the 
frequency of the first modified mode of the porous sheet. The lowest frequency boundary is 
naturally the solid-fluid decoupling frequency or at least the Biot frequency. 

Two absorption peaks close to unity were found for a H = 3.5 cm thick porous sheet, when 
a second circular cylinder of radius R^ = 5 mm is added to the configuration CI, figure [51 
The center of this cylinder is such that r\ = d = 2 cm and Q\ = n/3. The configuration C2 
was derived from a triangular lattice by reducing the radius of the upper cylinder to decrease 
the structure thickness. The first absorption peak stands around v\ R3 1850 Hz just below 
the first quarter-wavelength resonance frequency of the H = 3.5 cm thick porous sheet, 

(2) 

and the second stands around v\ ~ 4120 Hz. These two peaks correspond to minimum 
of |det (D) | and can therefore be explained by excitation of trapped modes, shifted in the 
complex plan. 

When the inclusions are inversely placed, i.e. the center of the first inclusion is 
(1 cm, 2.5 cm) and r^ 12 ) = d = 2 cm and 9 12 = — 7r/3, results are not identical and no 
quasi-total absorption peak is encoutered. This means that the configuration are not re- 
versible. 

The procedure was ran a second time with the addition of a third inclusion of radius 
R( 3 > = 2.5 mm to the configuration C2, r| = d = 2 cm and 9% = 27r/3. Three quasi- 
total absorption peaks were encoutered around v\ « 1500 Hz, v\ rs 3300 Hz, and v\ rs 
5000 Hz. This phenomenon was already encountered in 2 ^, where N trapped modes were 
found when iV cylinders were place across a wave tank. Nevertheless, this configuration 
imposes use of H = 5 cm thick plate and the absorption gain was considered unsignificant 
over the whole frequency range considered. 
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FIG. 5. Configuration C2 - Absorption coefficient of a H = 3.5 cm thick porous sheet of Firefiex 

backed by a rigid plate ( ) without inclusion embedded and ( — ) with a R 1 = 75 mm radius 

circular cylinder and a R 2 = 5 mm radius circular cylinder embedded with d = 2 cm. 



The addition of inclusions that imposes increasing the thickness of the structure is rapidly 
of non practical use because of the large absorption of the porous layer itself. 



Based on the fact that the frequency v\ decreases when x 2 increases, several attempts 
were followed to construct a porous sheet with a unit cell composed of varying x 2 central 
coordinate circular cylinder arranged in a kind of garland. For example, the absorption 
coefficient of a 8 circular cylinders per unit cell embedded in a 2 cm thick porous sheet 
was studied. The radius of the 8 cylinders was R^ = 75 mm and the projection of the 
center-to-center distance between two adjacent cylinders x^ was 2 cm. The x 2 were 

chosen such that x 2 — 1.1 cm, x 2 = 1.05 cm, x 2 — 1 cm, x 2 = 0.95 cm, x 2 = 0.9 cm, 
x 2 = 0.95 cm, x 2 — 1 cm, and x 2 = 1.05 cm. The absorption peak at v t was no more 
total and no particular increase of its width was noticed. This means that the periodicity 
has a large influence on the results and that the phenomenon can not be simply explained by 
trapped modes but rather by complex embedded Rayleigh-Bloch waves. A similar procedure 
was followed by decreasing the radius and the center-to-center distance, the cylinders being 
aligned, without particular effects on the absorption. 
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V. CONCLUSION 

The influence of embedding periodic circular inclusions on the absorption of a porous 
sheet attached on a rigid plate was studied theoretically and numerically. In addition to 
the absorption features related to excitation of modified mode of the plate and to Bragg 
interference, it is shown that the structure, in case of one array of cylinders embedded in 
a porous sheet, whose thickness and parameters, mainly the flow resistivity, are correctly 
chosen, possesses a quasi-total absorption peak below the quarter-wavelength resonance 
frequency. This particular feature enable the design of small dimension absorption packages 
and was explained by complex trapped mode excitation and by increase of the pressure 
gradient inside the layer. This quasi-total absorption peak was validated by use of Finite 
Element method, thus validating the described method and results. 

In case of more than one circular cylinder per spatial period, it was found that iV c 
quasi-total absorption peak can be obtained for particular arrangement along the porous 
thickness, i.e. close to triangular lattice. Nevertheless, this rapidly leads to a large thickness 
of the stucture and the embedding of the additional inclusions become useless. Garland 
arrangment were also tested without particular effect, or at least without as spectacular 
effect as the one already observed for one inclusion per spatial period. 

The method offers an alternative to mult i- layering and double porosity materials for the 
design of sound absoprtion packages. 

Appendix A: Modal analysis of the configuration 

The modes of the configuration without inclusions embedded (i.e. a rigid porous layer 
backed with a planar rigid wall), whose dispersion relation is 

D l = a [0]i cos (4 1]i #) - i« [1]i sin [k^ ]l H\ = , (Al) 

cannot be excited by a plane incident wave initially traveling in the air medium^. Effectively, 
figure E] depicts the real and the imaginary parts of the roots clAu) = u>/k*, n Jcu) of Eq. 
(IA1I) . as calculated for a H = 2 cm thick porous layer, whose acoustical characteristics are 
those used in section IIVI Under the rigid frame assumption and for frequencies higher than 
the Biot frequency (and lower than the diffusion limit), a porous material can be considered 
as a modified fluid, its associated dissipation being considered as a perturbation of a fluid. 
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For Eg. (I A I j) to be true without dissipation, k 2 should be purely imaginary while k 2 should 
be purely real. Under the previous assumptions, this implies that Re ( c*,, j should stand in 
[Re (cW) , c^], i.e. \k\\ should stand in [fc'°l; Re (fc' 1 ')]. Or for a plane incident wave initially 
propagating in the air medium \k\\ is always smaller than k^. It is also necessary to note 
that in the diffusion regime, i.e. for frequencies largely below the Biot Frequency, any mode 
exists. This fact constitutes the major difference when compared with a traditional fluid. 
Effectively, largely below the Biot frequency, fcW is purely imaginary. This implies that k 



[i] 



2 



is also purely imaginary whatever the value of k\ and that D l never vanishes. 

When inclusions are periodically embedded in the porous sheet, the dispersion relation of 
the modes of the configuration is det (I — V(S + Q)) = 0. The roots of the latter dispersion 
relation are difficult to determine because of the complex nature of the matrix I — V (S + Q). 
Here, we focuse on the case of only one grating, i.e. N c = 1 in order to emphasis the 
excitation of the modified mode of the backed layer (MMBL), 2 . Proceeding as in^, an 
iterative scheme can be employed to solved (ITTj) . The equation is re- written in the form 
(1 — VlM.ll)Bl = VlJ^l + Vl J2iez-^ Ll ^ 1 ^ ~ ^ Ll )- ^ ne iterative scheme reads as 

( £?f } = V L F L /{l-V L M LL ) 



< 



B {n+l} 



f V L Y,M Ll B} n \l-5Ll) + V L ^ , (A2) 

V zez / 

^ / (1 - V L M LL ) 

from which it becomes apparent that the solution B^ 1 , to any order of approximation, is 
expressed as a fraction, the denominator of which not depending on the order of approxi- 
mation can become small for certain couples (ki q ,uj), so as to make B^ 1 , and possibly the 
field large. 

When this happens, a natural mode of the configuration, comprising the inclusions and 
the plate, is excited, this taking the form of a resonance with respect to B^ , i.e., with 
respect to a plane wave component of the field in the plate relative to the inclusions. As 
B^ is related to f p , g p and R p , the structural resonance manifests itself for the same {k\ qi w), 
in the fields of the plate and in the air. 

The approximate dispersion relation 

V L = l-V L (s + J2QLL q ) =0, (A3) 
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is the sum of a term linked to the grating embodied in VlSq with a term linked to the plate 
embodied in Vl J2 q& zQLLq, whose expressions are give in (jHJ). 

This can be interpreted as a perturbation of the dispersion relation of the gratings by 
the presence of the plate. The zeroth order lattice sum can be re- written^ in the form 
Ylqzz^/dk^q (additional constants are neglected). Introducing this expression into ( 1A3j) 
gives 



ggZ ah '2q 1J q 



with 



M Lq = a[ 1] cos (k$H\ - ia [ q 0] sin fe ] #) + a[ 1] cos fc 1 (2x^ - if ) - 2L^ 

+ iaf sin (kf q hxP - H) - 2L9 q ) . (A5) 

It is then convenient, for the clarity of the explanations, to consider i) M^ to be a 
non-dissipative medium (a perfect fluid) and ii) the low frequency approximation of Vl, 
valid when A;' 1 ^ 1 ) « 1. The latter hypothesis ensures that the Vl reduces to VJ ~ 
(-l) i 7r(A;Wi?( 1 )) 2 /i4 + 0((A;W J RW) 2 ), I = -1,0,1. Equation (TM]) then reduces to 



■lYikWltt 1 



))2 



VimI-Y V ~ } ; ^ = , / = -1, 0, 1 . (A6) 

^ 2idk^D q /M q 

By referring to Cutler mode 24 , but also to the modal analysis carried out ini 4 -, the latter 
dispersion relation is satisfied (in the non-dissipative case) when the denominator of Eq. (1A6j) 
is purely imaginary and vanishes. These conditions are achieved when \ki q \ G [k'- ', Re (fc' 1 ')] 
and when either D q = or a q =0 (i.e. k\ = 0), which respectively corresponds to modified 
modes of the backed-layer (MMBL) and to modes of the grating (MG). Both of them are 
determined by the intersection of c\ q = iOJk\ q respectively with cIAuj) as calculated for the 
backed-layer and with Re (cW). The MMBL are pointed out by the dots on figure [6j The 
associated attenuation of each mode can then be determined by the values of Im ( tfs ) and 
Im (cW) at the frequencies at which the modes are excited. The attenuation associated with 
MG is also higher than the one associated with MMBL for all frequencies. Moreover, MG 
corresponds to the highest boundary of \ki q \ for Eq. (1A6l) to be true. This implies that MG 
should be difficult to excit. The latter type of mode can only be poorly excited by a plane 
incident wave, particularly at low frequencies. 
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FIG. 6. Real and imaginary part of the root of the dispersion relation in absence of inclusions C/ %, 
n = 1,2. Real part of the modified mode of the backed layer c? -,, n = 1, 2, q = 1, 2, for d = 2 cm 
are pointed out by dot. 
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